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Abstract 

It is shown that the momentum density of free electromagnetic field splits into two parts. One has 
no contribution to the net momentum due to the transversality condition. The other yields all the 
momentum. The angular momentum that is associated with the former part is spin, and the angular 
momentum that is associated with the latter part is orbital angular momentum. Expressions for 
the spin and orbital angular momentum are given in terms of the electric vector in reciprocal 
space. The spin and orbital angular momentum defined this way are used to investigate the 
angular momentum of nonparaxial beams that are described in a recently published paper [Phys. 
Rev. A 78, 063831 (2008)]. It is found that the orbital angular momentum depends, apart from 
an ^-dependent term, on two global quantities, the polarization represented by a generalized Jones 
vector and a new characteristic represented by a unit vector I, though the spin depends only on the 
polarization. The polarization dependence of orbital angular momentum through the impact of I is 
obtained and discussed. Some applications of the result obtained here are also made. The fact that 
the spin originates from the momentum density that has no contribution to the net momentum 
is used to show that there does not exist the paradox on the spin of circularly polarized plane 
wave. The polarization dependence of both spin and orbital angular momentum is shown to be the 
origin of conversion from the spin of a paraxial Laguerre-Gaussian beam into the orbital angular 
momentum of the focused beam through a high numerical aperture. 
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I. INTRODUCTION 



The orbital angular momentum (AM) of light did not draw much attention until 
1992 when Allen and his co-researchers 3j showed that a beam of Laguerre-Gaussian mode 
can carry both spin and orbital AM. They found that the spin is carried by the polarization 
a and the orbital AM is carried by the helical wave front represented by a phase factor 
exp(i/0), where I is an integer. Since then great progress has been made |4| in experiments. 



6|. The 



;ransfer of spin and orbital AM to microscopic 
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12 1 has been observed. 



The orbital AM has been measured 
particles BBS , [h]] and to liquid crystals 

Recently, experimentalists 3 showed that the spin and orbital AM of a non-paraxial 
beam play distinct roles in the interaction with microscopic birefringent particles trapped 
off the beam axis in optical tweezers. It was observed |l4| that the spin of light makes the 
particle rotate around its own axis and the orbital AM makes the particle rotate around the 



beam's axis. Furthermore, partial spin of a paraxial beam was observed [15| to be converted 
into orbital AM of a non-paraxial beam by a high numerical aperture. Those experimental 
results demonstrate that the spin and orbital AM of a non-paraxial beam are different in 
nature on one hand and are connected somehow to each other on the other. But up till now, 
there is no satisfactory theory to elucidate the difference and relation. The distinction that 



the spin is carried by the polarization and the orbital AM is carried by t 



re he 
3 



ical wave front 
, Q. It is not 



was drawn basically from the knowledge of a type of paraxial beams 
valid for non-paraxial beams [l8, Ji], 20]. With a specific non-paraxial beam, Barnett and 



Allen [18] found that "the seemingly natural separation of the angular momentum... is no 
longer possible" . The purpose of this paper is to advance a theory to explain the difference 
and relation between the spin and orbital AM of nonparaxial beams. 

To this end, we should first know how to represent a nonparaxial beam that as a whole 



is in a definite state of polarization. As mentioned before, Barnett and Allen 



put forward a nonparaxial solution. But that solution was shown [21] to fail to meet the 



181 ] once 



demand, 
developed 



brtunately, a theoretical representation that meets the demand was recently 



221 ] . The beam in this representation exhibits as a whole a definite polarization 



in the sense that all the plane waves that constitute the beam are described by the same 
normalized Jones vector. In other words, the normalized Jones vector in this representation 
is a global characteristic that plays the role of describing the polarization of the beam. This 
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Jones vector will be referred to as the generalized Jones vector. Apart from the global 
generalized Jones vector, a non-paraxial beam in this representation shows another global 
characteristic denoted by a unit vector. The global unit vector was applied {^J to explain 



the spin Hall effect of light 24|. In this paper, I will make use of this representation to 
show how the orbital AM depends on the polarization through the impact of the global unit 
vector. 

Secondly, we should also know how to define the spin and orbital AM of an electromag- 
netic field in free space. The total AM J(x ) of a free electromagnetic field with respect to 



the point x is defined as [25] 



J(x ) = J jd 3 x = 3(0) - x x J pd 3 x, (1) 

where j = (x — x ) x p is the AM density with respect to the same reference point, p = 
Eq/j,qSxT-[ is the momentum density defined in terms of the electric vector S and the magnetic 
vector TC, and 

J(0) = / x x pd 3 x (2) 



is the AM with respect to the origin. The separation of total AM into spin and orbital AM 



was discussed before [25|, 
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by performing the integration in Eq. ([Tj) by parts and 
neglecting a surface integral at infinity. In this paper, I will put forward a rigorous approach 
to the separation of total AM into spin and orbital parts by examining the property of 
momentum density. This approach allows us to apply the obtained result to plane waves. 

The paper is arranged as follows. In Section [Til h is found from the transversality 
condition that the momentum density of an electromagnetic field in free space splits into 
two parts. One part does not have any contribution to the net momentum; the other part 
produces all the momentum. The AM that originates from the former part does not depend 
on the choice of the reference point and is the spin. The AM that originates from the latter 
part is in general dependent on the choice of the reference point and is the orbital AM. In 
Section IHIt the integral expressions for the spin and orbital AM obtained in Section [TTJ are 
used to investigate the AM properties of nonparaxial beams described by the aforementioned 
representation. Since the light beam is assumed to be monochromatic, both the integrals 
of spin and orbital AM are infinite. In order to deal with the infinity, the technique of 
^-function normalization is used. As expected, the spin AM is found to be dependent on 
the polarization. But what is surprising is that the orbital AM is also dependent on the 



polarization. It is shown how the orbital AM depends on the polarization through the impact 
of the global unit vector. Two different problems are discussed in Section HVl by making use 
of the obtained results. Section [V] concludes the paper with further remarks. 



II. SEPARATION OF THE TOTAL AM INTO SPIN AND ORBITAL AM 

Consider an arbitrary electromagnetic field in free space. Its electric vector in real space 
can be expressed as an integral over the plane-wave spectrum, 

£(x,t) = l|^p J E(k)exp[*(k-x-u;i)]d 3 /c + c.c.j, (3) 

where k is the wave vector and E(k) is the electric vector in reciprocal space. The magnetic 
vector of the beam is derived from Eq. ([3]) and Maxwell's equations to be 



Integral expression (|3j) or (jlj) leads to the following transformations [25], 

u(-k) = -w(k), E(-k) = E*(k). (5) 

With the help of Eqs. (j3J) and (jlj) and vector algebra ax (b x c) = (a ■ c)b — (a ■ b)c, the 
momentum density splits into two parts, 

p = Eo/JoS x H = pi + p 2 , (6) 



where 



H 4(2tt) 3 J u 



A{2Ti) A J LO 

__EiL_ f ^^E*e^'-^ x e-^'-^d 3 k'd 3 k + c.c, (8) 
v ) J ^ 

E = E(k), E' = E(k'), uj = co(k), and u' = u(\s!). Based on the transversality condition 
k • E = 0, it is readily proven by use of transformations (J5J) that p 2 does not have any 
contribution to the net momentum, 

p 2 d 3 x = 0. (9) 



This tells us a fact that all the momentum P comes only from pi, 



E* E 

p = p x = / Vl d A x = e I - kd 3 k, (10) 



which is independent of time. 

Accordingly, the total AM also splits into two parts, 



J(x ) = y (x - x ) x pd 3 x = S(x ) + L(xq). 



Because of property (Q, the first part S that originates from momentum density P2 is 
independent of the choice of the reference point, 

S(x ) = S(0) = J x x p 2 d 3 x. (11) 

In other words, the fact that S is independent of the choice of the reference point roots in 
an intrinsic property of the electromagnetic field, the transversality condition. It is thus 
reasonable to regard this intrinsic AM as the spin. The second part that originates from 
momentum density pi is in general dependent on the choice of the reference point, 

L(x ) = L(0)-x xP 1 , (12) 

where 



L(0) = / x x Pl d 3 x. (13) 

It is plausible to regard this part as the orbital AM. Substituting Eq. (jSJ) into Eq. ffTTl) . one 
obtains by straightforward calculations 

S = f £^E* x Ed 3 k. (14) 
J ioj 

The momentum density P2 leads to the spin AM, though it does not produce any momentum. 
Such an astonishing fact means that there is no paradox on the spin AM of circularly 
polarized plane waves. This will be discussed in Section [TV] Substituting Eq. (J7J) into Eq. 
( TLB"]) , one has 

L(0) = / — E f (k x V k )Ed 3 fc, (15) 

J IU) 

where Vk is the gradient operator with respect to k, and the superscript f stands for the 



conjugate transpose 



291 ] . For the readers' convenience, the details to calculate Eqs. ([T 



and (|15p are summarized in Appendix. It is very interesting to note that the spin (Till) and 



orbital AM (ITBT) obtained this way look very like their quantum-mechanical counterparts 



25]. 



At last, let us give here for later convenience the total energy of the beam in terms of the 
plane-wave spectrum, 



(16) 



III. AM PROPERTIES OF NON-PARAXIAL BEAMS 



The AM o 
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a propagating beam in the z-direction is commonly considered in the literature 



19] to be equivalent to the line density, that is to say, to the AM per unit 



length in the z-direction. In order to avoid any possible ambiguity that may arise from the 
AM density flQ, I do not use th s notion here. In fact, we have given in Eqs. ( TH1) and 
(fl5l) the expressions for the spin and orbital AM themselves with respect to the origin. In 
this section, we will use those expressions to investigate the AM properties of nonparaxial 
beams. To do this, let us now convert the representation form of nonparaxial beams that 



was advanced in Ref. 



221 ] into a form that is suitable for present purpose. 



A. Description of non-paraxial beams: introduction to a new global unit vector 

The electric vector 8 of a nonparaxial beam in real space is given by Eq. . The electric 



22 



vector E in reciprocal space is factorized into three factors 

E = maf, (17) 

where 

m = (u v) (18) 
is the mapping matrix, a = ( ] is the generalized Jones vector that is assumed to be 

w 

independent of the wave vector and to satisfy the normalization condition a^a = 1, and / is 
the electric scalar in reciprocal space. The unit column vectors u and v of m represent the 
two mutually orthogonal states of linear polarization and are defined in terms of the local 
wave vector k and a global unit vector I as follows, 

k k x I , . 
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which lead to an important normalization property of the mapping matrix, 



m m = 1, 



(20) 



where the superscript T denotes the transpose. Unit vector I can be specified by its polar 
angle and azimuthal angle $. For the sake of simplicity, let us assume I to lie in the plane 
zox, that is to say $ = 0. In this case, we have 



and the mapping matrix 



m = 



fc|k x l| 



1(0) = e x sin + e z cos 



(ky + k 2 z ) sin — k z k x cos 



kk y cos 



\ 



(21) 



/ 



—k y (k z cos + k x sin 0) k(k z sin — k x cos 0) 
^ (k 2 + ky) cos — k z k x sin — kk y smQ 

where |k x I| = [k 2 — (/c x sin0 + k z cos 0) 2 ] 1 / 2 . Due to the symmetry relation 1(0 + n) = 
—1(0), it is postulated throughout this paper that 

7T 



101 < 



(22) 



A monochromatic beam has a definite wavenumber. It is convenient to use spherical 
polar coordinates to express the electric scalar as 

5(k - k') j 



f 



k 2 



where < $ < n and < <p < 2-k. Since f{d, <p) is a periodic function of <p with period 27r, 
a physically allowed function has the following Fourier expansion: 

oo 
i=— oo 

In this paper, we consider only one term of the expansion and rewrite the electric scalar as 
follows, 

(23) 



/ = ^ r ^/,Wexp(i^) J 



A; 2 

where the angular-spectrum function fiiti) is assumed to be square integrable. In order to 
use the technique of ^-normalization, the complex conjugate of / is replaced with 

S(k - k") 



r = 



k 2 



-/,*Wexp(-<V). 



(24) 



For a beam that propagates in the z-direction, its angular-spectrum function satisfies 

/,(#) = for | < i? < 7T. (25) 

Furthermore, if the beam is well-collimated and thus can be paraxially approximated, \fi{$)\ 
is sharply peaked at $ = 0. The half width A$ of is the divergence angle of the beam. 

So obtained E guarantees that the field vectors £ and TC in Eqs. (|3J) and (@J satisfy 
Maxwell's equations. Now that unit real vectors u and v are orthogonal to each other, the 
a that is independent of the wave vector acts as a global characteristic to describe the inner 
degree of freedom of the beam, the state of polarization. We thus have two independent 
global quantities, I and a, to describe a beam. It should be pointed out that a physically 
allowed beam may be a linear superposition of a series of so described beam. They each 
have their own I and a. The beam that we will consider in this paper is assumed to have 
definite I as well as a. In the following, we will pay much attention to the effect of these 
two global characteristics on the orbital AM. Only the AM with respect to the origin will 
be considered. 



B. Orbital AM is dependent on I as well as a 

The longitudinal component of orbital AM with respect to the origin can be turned from 
Eq. (USD into 

L z = [ —EU~i^-)Ek 2 smMkdMip (26) 
J u dtp 

in spherical polar coordinates. Hereafter the symbol for the origin will be omitted for the 
sake of simplicity. By making use of Eq. ( TT7I) , one has 

E, Hf)=*WH§^/7 + rH|). (27) 

When property (|2"U1) is taken into account, straightforward calculations yield 

i . . dm <7 3 cos $ — cos 



m'(—i-^—) = — <5"3Cos'# + 



d(p 2 1 — cos G cos d — sin sin •& cos tp 

+ <h cos fl + cos 9 

2 1 + cos cos t? + sin sin i? cos cp ' 
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-i 

where 03 = | ] is the Pauli matrix. Substituting Eq. (1281) into Eq. ( 1271) and noticing 

1 

Eq. (1231) . one obtains 

E'(— Z — ) = (l — (TCOSV)J J 1 



<9yj 2 1 — cos O cos $ — sin G sin i? cos <y2 

a (cos?? + cos Q)f*f 

2 1 + cos cos d + sin sin d cos </? 

Substituting it into Eq. fl26|) and considering Eqs. (1231) and ()24p . one finds after performing 
the integration with respect to variables k and p 

L z = ^S(k-k') f \m\ 2 ^dd+ 2 -^5{k-k') 

f n { 1 / cos t? + cos cos?? — cos 6 \ A , „ , „ s „ „ , „ . . 

x/ S- i n 7^ + 1 s -cost? >\M#) 2 sin -&d&. 29 

7 \2 V|cos?? + cos0| |cos??-cos0|y j lJlK n v ; 

In obtaining Eq. (|29|) . the following integral formula is used: 

<lj ' * "a| 2 < 1). (30) 



o 1 + a cos x y/T 



a 



2 



Substituting Eq. (ITT)) into Eq. (IT6l) and considering Eqs. (|23|) . (1241 . and (1201) . one has for 
the total energy of the beam 

W= ^S(k- k') \ft(#)\ 2 sin $d#. (31) 
k Jo 

It is clear that the orbital AM per unit energy is 

Lz = l_ a g Q ( gjSt| + | C c o S stc°o S s g| ) - «**} ^sinM 

Next let us calculate the transverse component of orbital AM. The x-component is rewrit- 
ten from Eq. (1151) to be 

L x = - f —E^lUt^-) - k z (i^-)]Ek 2 smtidkdtfdip. (33) 
J uj y dk z dk y 

According to Eq. (fTTj) . one has 

E%(* A) _ M , JL )]E = tf[k y m T (i^) - k z m T A]af*f 

z y z y 

+ rMi^)-k,(i-^u (34) 
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When property (1201) is taken into account, straightforward calculations yield 

T dm rp dm A a 3 (cos ^ — cos 9) cot 6 

k y m \ % -^r)- k z m ( l Tn~) = ^3 sm ^ cos v? + — . . — 

c/A; z oky 2 1 — cos B cos v — sm B sin v cos 9? 

03 (cos-!? + cos 6) cot B 

2 1 + cos 6 cos $ + sin B sin $ cos 97 

Substituting Eqs. ( 1341) and ( l35i) into Eq. ( 1331) and considering the rotation symmetry of / 
in Eq. fl23|) . one obtains after performing the integration with respect to variables k and ip, 

TlEnCr , ^ r ( COS 1? + COS B COS $ — COS B \ . „ . „, „ T „ , , 

L x = —r£-8(k-k' cotB / 1 — - + - /, tf 2 sin». 36 

/ru; J \| cosw + cosB| | cos 17 — cos B | / 

In obtaining Eq. fl36l) . formula fl30l) is used. The x-component of orbital AM per unit energy 
is thus 

~ f n 1 ( costf+cosB I costf-cose \ I fl„(f\ 12 ,Q-7,,Q 
L x _ CT COt B JO 2 ^|costf+cosO| + Icos^-coselj 



W u f£\f l (#)\ 2 ahi&d& 

Similar calculations give for the y-component of orbital AM per unit energy 



(37) 



Eqs. ( I32l . (1371) . and (I38I) are valid for any physically allowed angular-spectrum function 
Remembering that the unit vector I lies in the plane zox, they show that as a vector 
quantity, the orbital AM with respect to the origin is located in the plane formed by I and 
the propagation direction for the rotation-symmetry electric scalar ( 1231) . Apart from an 
/-dependent term in the longitudinal component, the orbital AM is closely dependent on the 
polarization a through the unit vector I. 

For a beam propagating in the z-direction, property ff25l) is satisfied. Considering our 
postulation (T22il . Eqs. (J32l) and (|37l) bocome 

L z _ I a Fj 2 {\ (l + gfgg| ) ~ cos*} |/,W| 2 sin^ 



W u u /* /2 |/,(#)|2sin0dtf 

= acot B Jo V2 I (l + gjggr) LM^sm^ 
^ w / " /2 |^(^)| 2 sin^ 



(40) 

/ I f, / /Ml - ^iii 

respectively. Eq. PU|) indicates that if I is neither perpendicular nor parallel to the prop- 
agation direction, the transverse component of orbital AM does not vanish. Let us discuss 
the following three cases. 
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This is the case in which I is perpendicular to the propagation direction. The beam 



described in this case is uniformly polarized 



22j in the paraxial approximation in the tradi- 



tional sense |32j. In this case, Eqs. (1391) and (jlDjl become 

Lm = 1 + n V2 (l-cosW^)| 2 sin^9 
W ~ uj uj / ff/2 |/,(i?)|2sini?dtf 

% = 0, 

respectively, indicating that the transverse component vanishes and the longitudinal com- 
ponent depends on the polarization. It should be noted that the vanishing transverse com- 
ponent here is just with respect to the origin. With respect to any reference point that is 
not on the beam axis (the 2-axis), the transverse component is by no means equal to zero 
as is shown by Eq. ffT2l) . Furthermore, by making use of paraxial approximation in which 
cosi!) in the integrand of the numerator can be approximated by unity, cos?? ~ 1, Eq. (I4ip 
reduces to 

— = (42) 
W uj K J 

Only under so special conditions, is the longitudinal component of orbital AM approximately 
independent of the polarization. Eq. (H2l) is exactly the result that was obtained from the 
consideration of paraxial Laguerre- Gaussian beams 

2. 6 = 

This is the case in which the unit vector I is parallel to the propagation direction. The 



beam described in this case is known as cylindrical vector beam 
and (HUj) become 

L z I cr/ " /2 |/^)| 2 costfsintfcM 
W uj uj J-/ 2 I 1 2 gin ^ 

— = 

W ' 



33 



34( | . In this case, Eqs. 



(43) 



respectively. The transverse component vanishes too. But it is seen from Eq. (l4~3]) that 
even in the paraxial approximation, the longitudinal component is not independent of the 
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polarization and is given by 

^ = l - - °-. (44) 
W uj uj K ' 



3. |e|»A0 

A well-collimated beam has a very narrow divergence angle Ai9. This situation allows us 
to consider such a case in whichl 0| 3> A$ is satisfied. The refracted beam that occurred 
in the spin Hall effect of light [24J was proven [2] to belong to this category. In this case, 
we have cosi? — cos© > in the region in which | /;($)! is appreciable. Eqs. (1391) and (flOl) 
are thus approximated as 

L z I af^O- -cos^)|/ ; (^)| 2 sin^ 
W ~ w uj j^ ,2 \f l (d)\ 2 sinMd 
L x a cot 



w ' u • < 45 > 

respectively. The longitudinal component is almost equal to that in the case of |0| = |. 
But the transverse component is not equal to zero. Eq. (|45|) expresses a simple polarization 
dependence through the unit vector I. 

C. Spin is dependent only on the polarization 

Substituting Eq. (fT7]) into Eq. ([H]) and taking Eqs. (123|) and (124")) into account, one gets 

S = IS^ " k ' ] j\\fiW^ d ^- (46) 
It shows that the transverse component of spin vanishes. The longitudinal component is 
given by 

S z =^p^5(k-k') ri/ z (tf)| 2 costfsin». (47) 
«™ Jo 



Clearly, the spin AM does not depend on the unit vector I. From Eqs. (1471) and (13ip . it 
follows that the longitudinal component of spin per unit energy is 

S; = a /;|/K^)| 2 cos^sin^ 

W w ^ |/z(^)| 2 sin^^ ' 1 j 

which is valid for any physically allowed angular-spectrum function f^fi). For a paraxial 
beam, cos ■& « 1 holds and Eq. (T4H!) reduces to 
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This is what was obtained from the consideration of paraxial Laguerre-Gaussian beams 



D. Total AM 

The total AM is the sum of spin and orbital AM. Since the transverse component of 
spin vanishes, we discuss here only the property of longitudinal component of the total AM. 
Combining Eqs. (1321) and (jUD together, one has 

t i n fnh ( i cos ^ cos % + r gs 4 zssa lr > l |/tWI 2 sinW 

J z L O JO 2 \| cos w+cos ©| | cos w— cos ©| y \J<>\ J\ 

W = u + u f*\f l (0)\ 2 smdd& ' ( ' jU) 

It is instructive to note that J z does consist of two parts. One depends only on an integer 
/, and the other depends only on a. But the former is not the orbital AM, and the latter 
is not the spin AM. Eq. (150|) is valid for any physically allowed function fi($). When Eq. 
( I25p is taken into account for a beam propagating in the z-direction, it becomes 

J z _ I | cr r /2 K 1 + ^^l)l^l 2sin ^ 

W uj uj / o ff/2 |/,(#)|2sin0d# 
which clearly shows the impact of the unit vector I. If = 0, Eq. (|5ip reduces to 

— = -, (52) 

which is independent of the polarization whether the beam is paraxial or not. If |0| = | on 
the other hand, one gets from Eq. (I5TT) 

which is also valid beyond the paraxial approximation. Though the total AM exhibits so 
simple dependence on / and a, the first term - is not the orbital AM and the second one ^ 
is not the spin AM, unless the paraxial approximation holds. It will be shown in the next 
section that the polarization dependence of L z for a nonparaxial beam of perpendicular I is 
the basis of conversion from spin to orbital AM by a high numerical aperture. 

In summary of this section, I have shown that the orbital AM is closely related to the unit 
vector I. It is due to the impact of I that the orbital AM is dependent on the polarization. If I 
is parallel to the propagation direction, both the spin and orbital AM have only longitudinal 
components. They are all polarization dependent whether the beam is paraxial or not. But 
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the total AM does not depend on the polarization. To the best of my knowledge, this is the 
first time to give the AM expression of cylindrical vector beams. If I is perpendicular to the 
propagation direction, the spin and orbital AM also have only longitudinal components. But 
in the paraxial approximation, the orbital AM is nearly independent of the polarization and 
is equal to - , and the spin AM is nearly equal to - . If I is neither parallel nor perpendicular 
to the propagation direction, the transverse component of orbital AM is not equal to zero. 
Comparison with the result of Ref. js| indicates that the unit vector I of Laguerre-Gaussian 
beams is perpendicular to the propagation direction. 



IV. APPLICATIONS 



In this section, I will apply the results obtained before to discuss two different problems. 
One is the so-called paradox on the spin of circularly polarized plane wave. It will be shown 
that such a paradox does not exist at all. The other is the conversion of partial spin of a 
paraxial beam to the orbital AM of the focused beam through a high numerical aperture. 
The conversion will be shown to root in the polarization dependence of both spin and orbital 
AM. 



A. There is no paradox on the spin of circularly polarized plane wave 



The so-called par adox on the spin of circularly polarized plane wave has been the subject 



of discussion ]2i 



35 



36] ever since Beth 37]] experimentally demonstrated that a circular! 



30 



38 



39 



40 



4011. 



polarized plane wave carries spin AM h and was still investigated recently 
It states that because the electric and magnetic vectors of a circularly polarized plane wave 
are perpendicular to the wave vector, its momentum density must be in the propagation 
direction. As a result, the AM component in the propagation direction must be zero 42 ] 
due to the cross product of the position vector with the momentum density. This is contrary 
to Beth's observation. 

As we have shown in Section [Til the spin of an electromagnetic field in free space does 
not come from the part of momentum density that produces the net momentum. Instead, 
it originates from the other part of momentum density that does not have contribution to 
the net momentum. From this point of view, it follows that there is no paradox on the spin 
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of circularly polarized plane wave. After all, what is produced from the momentum density 
in the propagation direction is the net momentum. In order to elucidate that the spin does 
not originate from this momentum density, let us make use of Eq. (1141) to calculate the AM 
of a plane wave. 

The electric vector of a plane wave in reciprocal space is given by 

E = ma/o<5 3 (k-k'), (54) 

where k' is the wave vector of the plane wave. If Eq. (154"1) is substituted directly into 
Eq. (Ti"4"|) . an infinity will occur. To deal with the infinity, we make use of the technique of 
^-normalization as before by replacing E* with 

E* = ma*f*5 3 {k-k"). (55) 

Substituting Eqs. fl54"|) and fl55|) into Eq. fTi"4"|) . one gets 

S = -e |/o| 2 ^ 3 (k-k'). (56) 

Similarly, substituting Eqs. fl54|) and fl55|) into Eq. ffl"6]) . one has for the total energy of the 
wave 

W = £o|/o| 2 5 3 (k-k'). (57) 
It follows that the spin per photon in the plane wave is 

S k 

— hw = ha-, (58) 

which is entirely along the direction of wave vector k. For circular polarizations a = ±1, 
the spin AM per photon is ±h, which is in perfect agreement with Beth's experimental 
observation. This indicates that when one talked about the paradox on the plane wave's 
spin, he/she did not realize the role that the momentum density in Eq. ([8]) plays in the 
AM. It is very interesting to note that we arrive at the quantum feature 25j of photon's 
spin by a purely classical approach, from which one might appreciate the nonlocal property 
of the photon. Since the spin comes from the momentum density that does not produce any 
momentum on one hand and is stored in the whole real space over which the plane wave 
spreads on the other, it might be probable that the concept of photon's spin density in real 
space is physically meaningless [3o|. 
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B. Conversion from spin to orbital AM by a high numerical aperture 



The incident beam in the AM conversion experiment 15| is LGq, a Laguerre-Gaussian 
beam. So its unit vector I is perpendicular to the propagation direction and its parameter 
/ is equal to one, 1 = 1. Before focusing, the spin and orbital AM per unit energy of the 
paraxial beam in the propagation direction are approximately - and -, respectively, as Eqs. 
( I49p and fj42|) show. After focusing, the spin per unit energy of the non-paraxial beam is 
obtained from Eq. (148]) to be 



a 



Io /2 IflW cos # sin Vdti 



v f* /2 \fl{#)\ 2 sin$d-d 

indicating that only a fraction of the incident spin remains in the focused beam, where fi($) 
now stands for the angular-spectrum function of the focused beam. If the rest of the incident 
spin 

"71-/2 




Jo l//W| 2 cos^sin^ 



OJ 



fo /2 \fi($)\ 2 smM<& 

is converted into the orbital AM [43|, the orbital AM of the focused beam should be 

/ " /2 |/^)| 2 costfsin^ 




/; /2 |/^)| 2 sin?W 

This is just the result predicted by Eq. (|4ip . We thus explain the conversion from the spin 
to the orbital AM on the basis that the orbital AM can be dependent on the polarization. 
If a = —1, the orbital AM per photon is less than h. On the other hand, if a = 1, the 
orbital AM per photon is larger than h. The authors of Ref. Q put forward their own 
theoretical explanation based on the analysis of the longitudinal component of the focused 
beam's electric vector. Because the longitudinal component of the electric vector is not able 
to represent the whole beam, they failed to show how the orbital AM of the focused beam 
depends on the polarization of the incident paraxial beam. 



V. CONCLUSIONS AND REMARKS 



In conclusion, I put forward a rigorous approach to the separation of the total AM into 
the spin and orbital AM. This approach is based on the analysis of the momentum density. 
It was shown that the momentum density can split into two parts. One part that does 
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not produce any momentum corresponds to the spin. The other part that produces all the 
momentum corresponds to the orbital AM. The spin defined this way was applied to show 
that there is no paradox about the spin of circularly polarized plane wave. Apart from 
the conclusion that the spin is dependent on the polarization, I further showed that the 
orbital AM is also dependent on the polarization. The polarization-dependent orbital AM 
was applied to explain the experiment jl5] that converted partial spin of the paraxial beam 
LGg into the orbital AM of the focused beam through a high numerical aperture. 

The unit vector I was shown to have evident impact on the orbital AM. In the first place, 
Eqs. ( |32|) . ( |37|) . and ( 138]) show that the orbital AM is located in the plane formed by I 
and the propagation direction. Secondly, Eqs. (T2"8"j) and (I3"5"l) show that the polarization- 
dependent term of orbital AM is determined by the direction of I. when I is parallel to 
the propagation direction, the orbital AM is always dependent on the polarization. When I 
is perpendicular to the propagation direction, the orbital AM is almost independent of the 
polarization in the paraxial approximation. These phenomena may imply that the orbital 
AM is most connected with the polarization, the inner degree of freedom, when I is parallel 
to the propagation direction and is least connected with the inner degree of freedom when 
I is perpendicular to the propagation direction. In a word, the impact of I on the orbital 
AM may offer further insights into the nature of the AM of light. 
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APPENDIX: DERIVATION OF EQS. (Q3D AND (1151) 

Let us first derive Eq. (jl~5]) . Substituting Eq. (J7|) into Eq. (fl~3]) . one has 



L(0) = Li + L 2 + c.c. 



(A.l) 



where 




(A.2) 



and 




(A.3) 
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Upon integrating Eq. (1A.2I) over the real space and noticing the following properties of 
Dirac's 5 function and its first-order derivative, 

1 r°° i r°° 

5(t) = — / exp(iu>t)du, 5'(t) = — / uj exp(iut)duj, (A. 4) 

27T J.oo 27T 



one obtains 

Li = — 

Ai 



J (k y e z - k z e y )^-^e~ i{u ' +uj)t 5'(k' x + k x )5(k' y + k y )5{k' z + k z )d 3 k'd 3 k 
+ £ -p [ (k z e x - k x e z )^-^e-^' + ^5{k' x + k x )8'(k' y + k y )5(k' z + k z )d 3 k'd 3 k 

Al J UJ y » 

+ g / (^e y - fc tf e ai )^5e- i ( ty+w >*<y(A4 + + + k z )d 3 k'd 3 k. 

It is changed by eliminating the 5 functions into 

U = £ f. j (k y e z - Key f 1 ^ ~ kz) ' E e -*K"04.-*».-*.)^(A4 + k x )dk x d 3 k 



+ % I (k z e x - k x e z ) " K ^ "^'" e-^^-^M^ + ^^A; 



4i 













E( 


kxi 




kg) 


•E 






u 






E( 




— h 

fry, 


K) 


E 



+ g J ( kx e y - k y e x ) E ( kx > ^ k ^r^ e -i[u>{- kx ,-k y K)Mt 5 ^ k ' z + k z )dk' z d 3 k. 



Noticing the following property of the derivative of the 5 function, 

,-t 2 



f(t)S'(t - t )dt = -/'(to), h < t < t 2 , (A.5) 



ti 



and taking transformation (jSj) into account, the above equation is reduced to 

L| = — / ^- ( E ■ ^ + z-^-E* • e] d 3 k 

Ai J u V 9k x SqUoU) J 

8W kyt , \ O, 

E • — — + i^^E* • E ) d k 



k y G z 




k z G y 




UJ 




k z G x 




k x G z 




uJ 




k x ^y 




k y &x 


UJ 



Ai J uj \ dk y Eo^quj 

4? / w V d/c 2 e /Jo^ 



^ / -E T (k x V k )E*d 3 k. 
A J uj 



By making the variable replacement k — > — k, it is changed into a familiar form, 

Li = I / ^Et(k x V k )Ed 3 fc. (A.6) 

Since operator — zVk is Hermitian, the hi in Eq. (1A.6I) is real. A similar calculation produces 
from Eq. fTA~3l) 

L 2 = L x . (A.7) 
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It is clear that substituting Eqs. flA.6jl and (IA.7D into Eq. (lA.ip will yield Eq. (1151) . 
Then we derive Eq. ( 1T41) . Substituting Eq. (jSJ) into Eq. (1TT1) . one has 



where 



and 



Si 



4(2tt) ; 



"4(2tt) j 



S = Si + S 2 + c.c, 



d 3 fc'(i 3 fc / d 3 x— -x x Ee i(k ' +k) - x e^ (tJ ' +tj)t ! 



d 3 k'd 3 k I d 3 x— -x x EV^'-k^e^'^*. 

a; 



(A.8) 



(A.9) 



(A.10) 



Upon integrating Eq. (1A.9|) over the real space and noticing Eq. (1A.4[) . one obtains 



Si 



l£ 



(EyQ Z E Z Gy 



+ — I (E z e x - E x e z 



Z£0 

4 

+ ^ / (E x e y - E y e x 



E'- 


k 


it) 


— 


E'- 


k 






E'- 


k 





e-i^-H^jfc/ + k x )S(k' y + k y )5{k' z + k z )d 3 k'd 3 k 
e-^'+^Sik'x + k x )S'(k' y + k y )5(k' z + k z )d 3 k'd 3 k 
e -i(a>'+a>)t 5 ^ + ^(jfc/ + ^ + k z )d 3 k'd 3 k. 



It is changed into, by eliminating the 5 functions and taking Eqs. (1A.5I) and (j5J) into account, 

<9E* 



s ' = f 



ie /" E y e z -E z ey ^ dE* ^ 3fc | ze /* £ z e x - g^e^ dE*^ 



4 



+ 



4 J " <% 2 

From the transversality condition k • E* = 0, we know that 

, dE* „ , dE* . , <9E* 
k . = k . = — _E* 7 k ■ 

dky 



-E* 



Si then reduces to 

Si = \ [ — E* x Ed 3 k, (A.ll) 
which is clearly real. Similarly, S 2 in Eq. ( 1A.10I) is found to be real and is equal to Si, 

S 2 = Si. (A.12) 



Substituting Eqs. fTATTT]) and (EH into Eq. (jOjl will yield Eq. (E 
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